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Abstract. Neutrino Deep Inelastic Scattering (DIS) on nuclei is an essential process to constrain the strange quark parton 
distribution functions (PDF) in the proton. The critical component on the way to using the neutrino DIS data in a proton 
PDF analysis is understanding the nuclear effects in parton distribution functions. We parametrize these effects by nuclear 
parton distribution functions (NPDF). Here we compare results from two analysis of NPDF both done at next- to-leading 
order in QCD. The first uses neutral current charged-lepton (i'^A) Deeply Inelastic Scattering (DIS) and Drell-Yan data for 
several nuclear targets and the second uses neutrino-nucleon DIS data. We compare the nuclear corrections factors (F^'°/F^) 
for the charged-lepton data with other results from the literature. In particular, we compare and contrast fits based upon the 
charged-lepton DIS data with those using neutrino-nucleon DIS data. 
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INTRODUCTION 

Parton distribution functions (PDFs) are an indispens- 
able part of any prediction involving hadrons in the ini- 
tial state. This is the reason why many groups perform 
and regularly update global analysis of PDFs for pro- 
tons [1, 2, 3, 4] and nuclei [5, 6, 7]. Although not of- 
ten emphasized, nuclear effects are present also in the 
proton PDFs analysis as a number of experimental data 
is taken on nuclear targets. Mostly though, the nuclear 
targets used in the proton analysis, are made of light nu- 
clei where nuclear effects are generally small. An impor- 
tant exception is the neutrino DIS data which is taken on 
heavy nuclei such as iron or lead and is sensitive to the 
strange quark content of the proton. A knowledge of the 
strange quark PDF has an influence on precise measure- 
ments at the LHC such as W- or Z-boson production. 

In order to make use of the neutrino DIS data to con- 
strain the strange quark PDF, we have to apply a nu- 
clear correction factor which can be obtained either from 
a specific model of nuclear interactions [8] or from an 
analysis of nuclear parton distribution functions (NPDF) 
based on experimental data. 

Here, we present a framework for a global analysis of 
nuclear PDFs at next-to-leading order in QCD closely 
related to the CTEQ framework for proton PDFs. We an- 
alyze and compare the nuclear correction factor obtained 
from the usual charged lepton DIS and Drell-Yan (DY) 
data to the one from the neutrino DIS data. 



NUCLEAR PDF 

The global NPDF framework, we use to analyze charged 
lepton DIS and DY data and neutrino DIS data, was 



introduced in [9]. The parameterizations of the parton 
distributions in bound protons at the input scale of Qq = 
1.3GeV 

xfk{x,QQ) = cox'i (1 -xy^e'^-''{l +e'*xy^ , (1) 
where k — Uy,dv,g,i{ + d,s,s and 

dix,Qo)/u{x,QQ) = cox'> (1 -x)"- + (1 +C3x)(l -xf* , 

(2) 

are a generalization of the parton parameterizations in 
free protons used in the CTEQ proton analysis [10]. To 
account for a variety of nuclear targets, the coefficients 
Cjt are generalized to functions of the nucleon number A 

Ck ^ Ck{A) = Ck.o + Ck.i ( 1 -A-'^'^) ,fe = {l,...,5}. 

(3) 

The proton PDF in this framework are obtained as a 
limit A — 1 and are held fixed at values obtained in 
the analysis [10]. From the input distributions, we can 
construct the PDFs for a general (A, Z) -nucleus 

/(^■^) e) = I Z;"/^ (X, Q) + X"/^ Q) , (4) 

where we relate the distributions of a bound neutron, 
f"^^ (x, Q), to those of a proton by isospin symmetry. 

We performed a global analysis of nuclear charged 
lepton DIS and DY data within this framework, deter- 
mining the A-dependence of the parameters c^^(A). In the 
analysis, we applied the same standard kinematic cuts 
Q > 2GeV and W > 3.5GeV as in [10] and obtain a fit 
with X^/dof of 0.946 to 708 data points with 32 free pa- 
rameters (for further details see [9]). 

The nuclear effects extracted in the form of NPDF are 
usually presented in the form of nuclear correction fac- 
tors. We focus on two nuclear correction factors related 




FIGURE 1. Nuclear correction factors i?^^(F2;ji:, Q^) (left) and (F2;x, Q^) (right) for global fits to charged lepton DIS and 
DY data (solid line) and to neutrino DIS cross-section data (dashed line) at the scale — 5GeV^. 



either to the DIS structure function F2 in the charged- 
current (CC) vA process 

or to the DIS structure function in the neutral-current 
(NC) l^A process 

[^^ + J'^ + ...]+4[m-^ + m'^ + ...] 

The superscript '0' stands for using the free nucleon 
PDFs ff[x,Q) and f^ix,Q) in Eq. (4) instead of the 

bound nucleon distributions fl'^^{x,Q) and f"^^{x,Q). 

In Fig. 1 (solid line), we show how the result of 
our global analysis of charged lepton data translates 
into these nuclear connection factors and how it com- 
pares to experimental data. As first observed in [11], the 
R'^^{F2;x,Q^) correction factor calculated using Eq. 5 
with parton densities from the fit to the charged lep- 
ton nuclear data, does not describe the NuTeV data 
well which raises the question if including neutrino DIS 
data in the global analysis corrects this behavior without 
spoiling the R^^'^{F2;x,Q^) correction factor which fits 
the charged lepton DIS and DY data well. 

NPDF FROM NEUTRINO DIS DATA 

To investigate the apparent discrepancy between the pre- 
dicted nuclear correction factor R^fj[F2',x,Q^) from the 
fit to charged lepton data and the neutrino charged cur- 
rent DIS data, we have set up a global analysis where 
we used exclusively the neutrino DIS cross-section data 
coming from NuTeV and Chorus experiments taken on 
iron and lead respectively. Here we applied the same 



TABLE 1. Summary table of a family of compro- 
mise fits. 



w 




xIa (/Pt) 


total ;t;^(/pt) 





638 (0.90) 




638 (0.90) 


1/7 


645 (0.91) 


4710(1.50) 


5355 (1.39) 


1/2 


680 (0.96) 


4405 (1.40) 


5085 (1.32) 


1 


736 (1.04) 


4277 (1.36) 


5014 (1.30) 


00 




4192 (1.33) 


4192(1.33) 



kinematic cuts as in the first analysis of the charged lep- 
ton data and we obtain a fit to 3134 neutrino DIS cross- 
section data points with j^/dof of 1.33 with 34 free pa- 
rameters (for further details see [12]). 

As was expected, the global fit to neutrino DIS data 
describes the data for the charged cuiTent nuclear cor- 
rection factor R^fj[F2\x^ Q^) well and does a poor job to 
describe the neutral current coiTection factor especially 
at low and intermediate Bjorken x. We see that using one 
or the other data sets produces different nuclear correc- 
tion factors R''i^^{F2;x,Q^) and R^f^{F2;x,Q^). A ques- 
tion arises if there are such nuclear correction factors 
which would be in agreement with both charged lepton 
and neutrino data, for example using a combined set of 
charged lepton and neutrino DIS data. Analyzing both 
data sets in a combined global analysis runs into the 
problem of imbalance of number of data points between 
the two data sets. This would automatically mean that the 
neutrino data would be favored just based on the amount 
of data. Therefore, we introduce an artificial parameter, 
the weight of the neutrino data set w, 

X'^ E Xf + L ""Xf , (7) 

/±A data vA data 

to interpolate between the two different global fits (w = 
results in the fit to charged lepton data only and w ~ °° 
stands symbolically for the fit only to neutrino data). 
Varying the weight w, we try to find a compromise fit 




X 



FIGURE 2. Nuclear correction factors i?^^ (^2 ; ^, 2^ ) (left) and {F2;x,Q^) (right) for compromise fits with different weights 
of the neutrino DIS data at the scale = SGeV^. 



which would describe both charged lepton and neutrino 
data well. We list the resulting for the compromise 
fits with weights w = 0, 1 /7, 1 /2, l,oo in Tab. 1 and we 
show the nuclear correction factors R^^'^{F2;x,Q^) and 
R^^[F2\x,Q^) for the same family of compromise fits in 
Fig. 2. 

Indeed we see in Fig. 2 that the fits with w = 
1/7,1/2,1 interpolate well between the two extreme 
cases w = and w = °°. In order to decide on how well 
the compromise fits describe the data we use the 
goodness-of-fit criterion introduced and used in [13, 2]. 
We consider a fit a good compromise if its for both 
data subsets, the charged lepton DIS and DY data and 
the neutrino DIS data, is within 90% confidence level of 
the fits to only charged lepton or neutrino data. 

We define the 90% percentile ^90 used to define the 
90% confidence level, by 

/ P{x\N)dx^ =0.90, (8) 
Jo 

where is the number of degrees of freedom and 
P{X^,N) is the probability distribution 

P(7 N) = (9) 

' ' 2^/2r(A^/2) ■ ^ ' 

We can assign a 90% confidence level error band to the 
X^ of the fits to the charged lepton DIS and DY data and 
to the neutrino DIS data 

;t,l^ = 638 +45.6, = 4192+ 138. (10) 

Comparing the results of the compromise fits with differ- 
ent weights, listed in Tab. 1, we see that none of the com- 
promise fits are compatible with both 90% confidence 
level limits given in Eq.lO. As detailed in [12], not even 
relaxing the condition to compare against the 99% confi- 
dence level limit helps to finding a suitable compromise 
fit. Moreover, we show in [12] that the effect is related to 
the precise neutrino DIS data from NuTeV. 



CONCLUSION 

After performing a thorough global NPDF analysis of 
the combined charged lepton and neutrino data, we find 
that there is no good compromise description of both the 
data sets simultaneously. The differences are most pro- 
nounced in the low and intermediate x regions where the 
neutrino DIS (NuTeV) do not show a strong shadowing 
effect as the charged lepton data do. The inability to de- 
scribe all data by one consistent framework indicates the 
existence of non-universal nuclear effects or unexpect- 
edly large higher-twist effects, here 
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